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I.  INTRODUCTION 


The  thrust  of  a  rocket  is  produced  by  the  high-speed  discharge  of 
gases  to  the  rear  of  the  rocket.  In  effect,  the  rearward  linear  momen¬ 
tum  of  the  rocket  gases  is  balanced  by  an  increased  forward  momentum 
of  the  rocket.  In  the  case  of  a  rocket  with  pitching  and  yawing 
motion,  these  gases  also  possess  angular  momentum  in  pitch  and  yaw  and 
a  careful  dynamics  analysis  shows  a  damping  of  the  pitching  and  yawing 
motion  associated  with  this  gaseous  angular  momentum.  This  pitch  jet 
damping  is  usually  fairly  small.  For  Army  rocket -assisted  projectiles 
(RAPs),  it  is  completely  negligible. 

A  spinning  RAP  ejects  gases  with  spin  angular  momentum  and  the 
effect  of  this  lost  angular  momentum  should  be  considered.  In  this 
report  an  expression  for  the  spin  jet  damping  of  a  rocket  is  derived 
and  the  size  of  this  effect  is  computed  for  two  Army  RAPs:  the  155mm 
M549  and  the  8-inch  M650 .  The  predicted  change  in  spin  will  then  be 
compared  with  actual  flight  measurements  made  by  sunsondes  with  onboard 
telemetry  transmitters1"2. 


II.  THEORY 

We  will  use  the  usual  nonrolling  aeroballistic  x,  v  z  axes  with 
the  x-axis  along  the  rocket's  axis  of  symmetry  and  the  z-axis 
selected  to  be  initially  downward -po int ing .  If  the  rocket  is  flying  at 
zero  angle  of  attack,  the  only  aerodynamic  force  acting  on  the  rocket 
is  the  drag  force  ?  .  For  small  angles  of  attack  and  sideslip,  the  x- 

component  of  this  force  is  F^  =  |  .  The  x-component  of  the  linear 

momentum  of  the  rocket  body  plus  propellant  at  time  t  will  be  denoted  by 
A  (t).  The  derivative  of  A^  will  then  be  equal  to -F^  plus  the  x- 

component  of  the  gravity  force. 


1.  William  H.  Mermagen,  Measurements  of  the  Dynamic  Behavior  of 
Projectiles  over  Long  Flight  Paths , 99  Journal  of  Spacecraft  and 
Pockets  8,  April  1971*  pp.  380-385 .  (See  also  Ballistic  Research 
Laboratories  Memorandum  Report  No.  2079 ,  November  1970 ,  AD  717002.) 

2.  Charles  H.  Murphy,  nEffect  of  Large  High-Frequency  Angular  Motion 
of  a  Shell  on  the  Analysis  of  Its  Yawsonde  Records,  rr  Ballistic 
Research  Laboratories  Memorandum  Report  No.  2581,  February  1976, 

AD  B009421L. 
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(1) 


It  is  very  important  to  note  that  the  derivative  in  Equation  (1) 
is  the  rate  of  change  of  linear  momentum  for  the  constant-mass  system 
of  rocket  plus  propellant  and  must  be  calculated  with  care.  At  time 
tQ  +  At  the  rocket  has  a  mass  m  +  Am  (where  Am  is  negative)  and  an  x- 

component  of  velocity,  u  +  Au.  A  mass  of  propellant  gas,  -  Am  ,  has 
separated  from  the  rocket  and  is  traveling  away  from  the  rocket  in  the 
exhaust  jet  with  a  velocity  x-component  Uj  relative  to  the  rocket. 

The  new  linear  momentum  is,  therefore. 


A  (t  )  +  AA  =  (m  +  Am)  (u  +  Au)  +  (-  Am)  (u  +  u  )  (2) 

X  0  X  J 


where 


A  ft  )  =  m  u 
x  0 


Equation  (2)  is  divided  by  At  and  the  limit  as  At  +  0  is  taken  to 
obtain  .  This  result  can  be  substituted  in  Equation  (1)  to  yield 


m  u  =  Ft  -  Fd  +  m  gx  (3) 


where 

F  =  m  u  is  the  x-component  of  the  thrust 

^  ^  of  the  rocket. 

This  simple  derivation  of  rocket  thrust  has  been  done  so  that  we 
can  see  the  proper  way  to  estimate  the  jet  effect  on  the  angular  motion 
of  the  rocket.  If  the  components  of  the  angular  momentum  of  the  rocket 
are  (H^,  H- ,  H-) ,  the  angular  motion  must  satisfy  the  following  vector 

differential  equation : 


H-,  H-)  +  (0,  q,  r)  X  fHx,  H~,  H~) 


CMX,  M~ ,  M.) 


(4) 
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We  will  assume  that  the  rocket  grain  burns  symmetrically  so  that  the  x- 
axis  remains  an  axis  of  mass  symmetry  and  the  transverse  moments  of  in¬ 
ertia  remain  equal.  The  angular  momentum  of  the  solid  part  of  the 
rocket  is  then 

Ss  =  (Ix  rt  rt  (5) 

and  at  time  t 

0 

ft  =  (H  ,  H~,  H~) 
x  y  z 

=  ax  p  +  hx,  It  q  +  h~ ,  It  r  +  h~)  (6) 

where  the  h.'s  are  the  angular  momentum  components  of  the  gas  in  the 
rocket  motor. 

At  time  t  +  At,  the  angular  momentum  of  the  solid  part  has  slight¬ 
ly  changed  and  some  propellant  gas  of  mass  -  Am  has  been  emitted  from 
the  nozzle,  but  the  small  change  in  the  angular  momentum  of  the  gas  in 
the  motor  will  be  neglected. 

The  angular  momentum  lij  of  the  emitted  gas  can  be  expressed  as  a 
sum  of  volume  elements: 

i\j  x  V  PJ  4(Vol)  t7) 

where  it  ,  it  and  pT  are,  respectively,  the  position,  velocity  and 
J  J  J 

density  of  the  jet  of  propellant  gas  in  the  particular  volume  element: 


it  =  (x,  R  cos0,  R  sin0) 

J 

(8) 

■  (uj>  fj>  "j1 

(9) 

A(Vol)  =  R  AR  A0  Ax 

(10) 
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and  where  R  and  0  are  polar  coordinates.  In  the  time  interval  At,  the 

emitted  gas  forms  a  disk  of  thickness  Ax  =  u  At  and  radius  equal  to 

j 

that  of  the  rocket  nozzle  exit,  R  ,  located  a  distance  X  from  the 
.  •  -i  n  n 

projectile  c.m. 

hj  =  Pj  Ax  jRn  |  * (R (Wj  cose  -  Vj  sine)  , 

o  0 

R  uT  sin0  -  X  wT  ,  (11) 

J  n  J 

X  vT  -  R  uT  cos0  ^R  dR  d0 
n  J  J  / 


The  transverse  motion  of  the  exhaust  gases  is  assumed  to  be  the  sum  of 
the  pitching  motion  of  the  nozzle  and  a  rotationally  symmetric  circum¬ 
ferential  motion  due  to  projectile  spin. 


(12) 

(13) 

where 

R  =  R/R 

n 

f(0)  =  0,  f(l)  =  1  . 


vT  =  r  X  -  pR  f  (R)  sinO 
J  n  r  n  v  J 

wT  =  -  q  X  +  p  R  f (R)  cose 
.T  n  n  r  n 


Equations  (12-13)  allow  the  angular  momentum  of  the  gas  emitted  in  the 
time  interval  At  to  be  written  in  a  very  simple  form. 


hj  =  (Rj  F  p,  X2  q,  X2  ^)  (-  Am) 


(14) 
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where 


,1 


F  =  2 


R2  f(R)  dR 


Equation  (14)  now  allows  us  to  compute  the  change  in  angular 
momentum  of  the  system. 


H  +  AH  =  (I  +  AI  ) (p  +  Ap)  +  h  -  (R2  F  Am)p  (15) 

x  x  x  xrr  x  n 


H~  +  AH~  =  (I  +  AI  ) (q  +  Aq)  +  h~  -  (X2  Am)q 
y  y  t  t  y  n 


(16) 


H~  +  AH~  =  (I  +  AI  ) (r  +  Ar)  +  h~  -  (X2  Am)r  (17) 

z  z  l  t  z  n 


Equations  (15-17)  can  be  simplified  by  Equation  (6),  divided  by 
At,  and  the  limit  taken  as  At  -*■  0  to  yield  the  derivatives  of  the 
angular  momentum  components.  Then  Equation  (4)  becomes 


I 

x 


P 


q  h~ 
z 


r  h~  =  J  p  +  M 
y  p  F  x 


(18) 


It(q  +  i  r)  -  i  (Ix  p  +  hx)(q  +  i  r)  -  Jq(q  +  i  r)  +  (M~  +  i  M~) 

(19) 


where 


J  =  m  R2  F  -  I 
p  n  x 


J  =  m  X2  -  i 
q  n  t 
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The  small  terms  in  the  angular  momentum  of  the  motor  gas  are  usually 
neglected  so  that  the  only  effects  of  the  propellant  gases  on  the 
angular  motion  of  the  rocket  are  the  time-varying  moments  of  inertia, 
the  spin  jet  damping,  J  p,  and  the  pitch  jet  damping,  J  (q  +  i  r)  . 


III.  DISCUSSION 


J 

q 


is  the  derivative  of  the  change  in  pitch  moment  of  inertia  of 


the  propellant  gases  and  is  negative.  Thus  this  term  damps  the  pitch¬ 
ing  motion,  producing  the  well-known  pitch  jet  damping  effect.  As  we 
shall  see,  for  Army  RAPs  the  pitch  jet  damping  moment  is  much  smaller 
than  the  aerodynamic  damping  moment  and  has  no  measurable  effect  on  the 
pitching  motion. 


If  we  assume  that  the  radial  variation  in  circumferential  velocity 
is  linear  like  that  of  a  spinning  rigid  body,  then  f(R)  =  R  and  F  =  1/2. 
For  this  case,  is  the  derivative  of  the  change  in  spin  moment  of 

inertia  of  the  propellant  gases.  Since  is  positive,  this  term 

represents  a  spin  jet  undamping.  Since  it  is  unreasonable  to  expect 
the  spin  angular  momentum  of  the  emitted  gases  to  drop  to  the  low  value 
associated  with  f(R)  =  R,  this  value  yields  a  maximum  value  of  and 

we  would  expect  the  actual  value  to  be  less. 


where 


(J  ) 
v  p  max 


0  <  c  <  1 


(J  ) 
p  max 


iii  (R2/2)  -  I 
^  n  '  x 


(20) 


The  aerodynamic  moment  contains  a  spin  damping  term  in  p  and  a 
pitch  damping  term  in  q  +  i  r.  In  order  to  compare  the  aerodynamic 
terms  with  the  jet  damping  terms  we  will  nondimensionalize  the  jet 
damping  moments  in  the  same  way  that  the  aerodynamic  moments  are  non- 
dimensional  ized  3  . 


3.  Charles  H.  Murphy,  "Free  Flight  Motion  of  Symmetric  Missiles ,  " 
Ballistic  Research  Laboratories  Report  No .  1216,  July  1963, 

AD  442757 . 
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(21) 


P 


J,  =(i)" vs  12  cj 


q 


(22) 


In  Table  1  the  average  values  of  I  ,  I  ,  m  and  other  appropriate 

X  L 


parameters  are  given  for  two  rocket-assisted  projectiles,  the  155mm 
M549  and  the  8-inch  M650.  According  to  this  table,  the  pitch  jet 
damping  coefficient  is  much  smaller  than  the  aerodynamic  pitch  damping 
coefficient  and  can  be  neglected,  but  the  spin  jet  damping  coefficient 
is  the  same  size  as  the  aerodynamic  spin  damping  coefficient.  We 
would  therefore  expect  the  spin  to  be  affected  by  this  term  during 
burning. 

Another  way  to  estimate  the  effect  of  the  spin  jet  damping  is  to 
consider  Equation  (18)  for  no  aerodynamic  moment  and  neglecting  h~  and 


y 


h~ . 


z 


(23) 


where 


y  can  be  computed  from  the  table  and  is  0.90  and  0.87  for  the  M549  and 
M650,  respectively.  Thus  the  percentage  change  in  spin  during  burning 
can  be  approximately  related  to  the  percentage  change  in  the  spin 
moment  of  inertia. 


Ap 

— =  -  e  y 
P 


(24) 


0 . 046e 


(M549) 


0. 033e 


(M650) 
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TABLE  1 .  PARAMETERS 

FOR  TWO  ROCKET -ASSISTED 

PROJECTILES 

155mm  M549 

8 -inch  M650 

burn-time  (s) 

2.68 

2.95 

m  (kg/s) 

-  1.12 

-  1.86 

Ix  (kg-m2/ s) 

-  .00285 

-  .00714 

It  (kg-m2/ s) 

-  .0379 

-  .0803 

Rn  C») 

.0222 

.0317 

Xn  W 

.326 

.365 

V  (m/s) 

481 

603 

c£ 

-  .0112 

-  .0113 

P 

CM  +  CM. 
q  ct 

-13.2 

-11.9 

(c  ) 

.019 

.012 

\  p/max 

CJ 

-  .61 

-  .34 

q 
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IV.  EXPERIMENTAL  RESULTS 


Recently  yawsonde  data  have  been  obtained4-5  for  the  spin  during 
burning  of  both  the  M549  and  the  M650.  Sample  spin  histories  for  each 
RAP  are  given  in  Figure  1.  The  spin  curve  before  burning  was  extended 
through  burning  and  the  percentage  change  in  spin  during  burning  was 
determined.  It  was  1.9%  for  the  M549  and  1.7%  for  the  M650.  These 
values  correspond  to  e  values  of  0.41  and  0.52,  respectively. 


V .  SUMMARY 

1.  The  derivation  of  rocket  thrust  and  pitch  jet  damping  has 
been  reviewed. 

2.  A  theoretical  model  for  spin  jet  damping  has  been  derived 
which  predicts  a  percentage  increase  in  spin  proportional  to  the  per¬ 
centage  decrease  in  the  spin  moment  of  inertia. 

3.  Experimental  results  show  approximately  2%  increase  in  spin, 
which  is  about  half  the  predicted  maximum  increase. 


4 .  Anders  S.  Platou ,  "Yawsonde  Flights  of  155mm  Non-Conioal  Boattail 
Progeotiles  and  the  155mm  M549  Progeotile  at  Tonopah  Test  Range — 
Ootoher  1977*1  USA  ARRADCOM  Ballistio  Researoh  Laboratory  Memorandum 
Report  No.  ARBRL-MR-02881 ,  November  1978,  AD  A065356. 

5.  Vural  Oskay ,  Wallaoe  H.  Clay  and  Martin  Klawa,  "8-Inoh  PXR  6263 
(Temperature  Instrumented  M650  RAP)  Tests  at  Yuma  Proving  Ground,  " 
USA  ARRADCOM  Ballistio  Researoh  Laboratory  Memorandum  Report  No. 
ARBRL-MR-03037 ,  August  1980. 
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Figure  1.  Measured  spin  histories  for  two  rocket-assisted  projectiles,  showing  the  percentage 
spin-up  during  burning. 
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A 

X 

x-component  of  the  linear  momentum  of  the  rocket 
body-plus-propellant  system 

CJ 

P 

j 

p _ 

(1/2)  p  V  S  £2 

CJ 

q 

J 

q _ 

(1/2)  p  VS  t2 

C£ 

P 

roll  damping  moment  coefficient 

CM  +  CM. 
q  a 

sum  of  the  damping  moments 

(1/2)  p  V2  S  £  (q2  +  r2)^ 

F 

fl 

2  R2  f  (R)  dR 

•’o 

fd 

l?Dl 

drag  force 

ft 

m  u  ,  x-component  of  the  rocket  thrust 

«j 

f(R) 

differentiable  function  of  R,  where  f(0)  =  0, 

f Cl)  =  i  . 

gx 

x-component  of  gravity  acceleration 

3 

angular  momentum  of  the  rocket 

fis 

angular  momentum  of  the  solid  part  of  the  rocket 

H  ,  H~,  H~ 
x  y  z 

aeroballistic  system  components  of  H 
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hJ 
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J 
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m  R^  F  -  i  =  (£/V)  •  spin  jet  damping  moment 

J 

q 

•  A  * 

m  -  I  =  (£/V)  •  pitch  jet  damping  moment 

£ 

reference  length 

M  ,  M-,  M~ 
x  y  z 

aeroballistic  components  of  the  aerodynamic  moment 

m 

mass  of  the  rocket  body-plus-propellant  system  at 
time  t 

0 

p»  q,  r 

rocket  spin,  pitch  and  yaw  rates  in  the  aero¬ 
ballistic  system 

R 

polar  distance  coordinate,  0  <  R  <  R^ 

A 

R 

position  vector  of  the  gas  in  the  exhaust  jet 

R 

n 

radius  of  the  rocket  nozzle  exit 

S 
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t 

time 

u 

x-component  of  the  velocity  of  the  rocket  body- 
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a 


V 

x,  y,  z 

X 

n 

Y 

Am 

e 

0 

P 


magnitude  of  the  velocity  vector 

velocity  of  the  gas  in  the  exhaust  jet  relative  to 
the  rocket 


nonrolling  aeroballistic  axes,  with  the  x-axis 
along  the  rocket's  axis  of  symmetry  and  the  z- 
axis  initially  downward 

distance  from  the  rocket  c.m.  to  the  rocket 
nozzle  exit 


the  change  in  mass  of  the  rocket  body-plus-pro¬ 
pellant  system  in  the  time  interval  At 

a  nondimensional  constant  in  the  range  0  to  1  . 

polar  angular  coordinate 

air  density 

density  of  the  gas  in  the  exhaust  jet 


V 

x,  y,  z 

X 

n 

Y 

Am 

E 

0 

P 


Subscript : 


( ) 


max 


value  at  f(R)  =  R,  F  =  1/2 
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1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


